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We discuss the possibility that fermions bind due to Higgs or pseudoscalar exchange. It is 
reasonable to believe on qualitative grounds that this can occur for fermions with a mass larger 
than 800-900 GeV. An exchange of a pseudoscalar boson leads in the nonrelativisitc limit to an 
unacceptable potential which behaves like 1/r 3 at the origin. We show that this singular behaviour 
is smeared out when relativistic effects are included. 



I. INTRODUCTION 



' In this note, we obtain an estimate for the mass of a fermion which is heavy enough to bind to its anti-particle via 
Higgs exchange. We then try to repeat the argument for pseudoscalar exchange, since pseudoscalars arise naturally in 
some extensions of the standard model. We will find eventually that the situation is similar to that in Higgs exchange, 
but the analysis is now much more complex as the exchange of a non-relativistic pseudoscalar seems to lead to a 
potential which is singular at the origin. We shall show that if the pseudoscalar exchange is treated in a relativistic 

| framework this singularity vanishes. 

We begin by considering Higgs exchange between the usual quarks taken to be of mass to. This leads to a weak 
attraction with coupling constant (see for example Reference [Q): 

/ rji \ / \ — 1/2 

t-h I g=V2(—)«l; v=(G F V2) = 245 GeV (1) 

Even for the t-quark with m t ~ 180 GeV, g ~ 1. (This does not necessarily invalidate a perturbative approach 
since the typical parameter of the perturbative expansion is g 2 /STT 2 <C 1). Even this coupling, however, is not strong 
i-S^ ' enough to create bound states of ti or tt. In principle, however, heavier fermions / for which the coupling g is greater 
than one can be found where bound states // are bound by Higgs exchange. 

At a qualitative level it is easy to estimate how heavy fermions should be for them to be bound by means of Higgs 
exchange. Obviously for non-relativistic particles, the kinetic energy p 2 /m (we use the reduced mass to/2 of the two 
particles with mass to) should be smaller than the potential energy which is 



\S\ = fle- m » r (2) 
47T r 

where ran is the Higgs mass. If the composite state has the radius a, then p 2 /to ~ 1 / ma, 2 and we get the condition: 



This equation can be rewritten in the form 



1 q 2 1 

< f — e~ mHa (3) 



g 2 m 4-Ke mHa 

> (4) 



niH mud 

Whatever the value of a the right hand side of Eq (0) is always larger than 47re. Therefore 



2 

^ > 47re = 34.3 (5) 
m H 
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We use Eq (|j) to estimate g. So 



to 3 
v 2 rriH 



> 17 (6) 



For to = m t = 180 GeV this leads to tuh < 5.7 GeV, which is unacceptable. However if we take to# = 400 GeV (for 
example), Eq. (|J) gives: 

m > 800 GeV (7) 

This would be a very heavy fermion, but fcrmions of mass in the TeV region arise in many extensions of the standard 
model. 

We will now try to extend this simple analysis to pseudoscalar exchange. Though in the SM there is no pseudoscalar 
Higgs they occur in many extensions of the SM, for example in supersymmetric extensions. We show in the next 
Section that the exchange of a pseudoscalar meson A leads in the nonrclativistic approximation to the potential 



2 



\ (<ti.(7 2 - 3(<Ti.n)((T2.n)) (1 + m A r) (8) 



■in 



(ai.n)(a 2 .n) 

r 

where a\,a2 are the Pauli matrices for the two fermions, n — rjr and tiia, va are the mass of the A-boson and the 
scale, analogous to v = 246 GeV in the standard model. The potential of Eq (j|) depends on the spins of the fermions 
and in certain spin states can lead to an attraction between fermions. In this case the analogue of the inequality (^) 
is 

< T^\ e ~ mAa ( 9 ) 
ma z 'vkv A or 

Clearly for a small enough radius a this condition can be satisfied. But this corresponds to the well-known fact that 
an attractive potential like 1 /r 3 which is more singular at the origin than the centrifugal term leads to the "collapse" 
of the particle to the origin and hence to a problem which is not well-defined Q . We shall demonstrate, however, that 
while the expression (||) holds for nonrelativistic fermions and leads to collapse, a relativistic treatment will lessen the 
singularity at r — so that there is no collapse to the centre. 

Suppose now that 

m A a < 1 (10) 

Then (H) means that 



The condition that the fermions should be nonrelativistic p 2 jm < m at p ~ 1/a together with eq. ( pi] ) can be written 
in the form: 



I TO 



< a <— (12) 



It is thus required that 

For va ~ v — 245 GeV this leads to 



to > v A-kva (13) 



to > 870 GeV (14) 

This restriction is not very different from (0) for the scalar case. Note that for a ~ (870 GeV) -1 the condition ( fTo| ) 
is probably satisfied since we can expect the pseudoscalar A to be lighter than the Higgs (the present limit is only 
ni a > 24.3 GeV) @. 

Note an interesting feature of bound states mediated by scalar and pseudoscalar exchange: unlike the situation with 
vector exchange, states of the same mass can appear for both // and // systems. This is a consequence of the fact 
that scalar and pseudoscalar interactions are even under charge conjugation. 
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II. PSEUDOSCALAR EXCHANGE (NON-RELATIVISTIC CASE) 



We now look in more detail at pseudoscalar exchange. Start with the interaction of the pseudo-scalar boson A with a 
fermion via the coupling i(2m/vA) yPlbW) A. The Feynman amplitude corresponding to the exchange of A can then 
be represented in the form 

4r7i 2 1 
M = —3- ((p 2 *aq(p 2 ) ((p*aqipi) 2 „ , (15) 

where q = — p\ — p 2 — j? 2 IS the three-momentum transfer and we have taken the non-relativistic limit so that the 
four-component spinors ip are replaced by the two-component spinors (p. 



The Born amplitude A (q) is related to M (q) by 



A = 



1 



-M 



8nW 

where W is the total centre of mass energy, and the potential P (r) is the Fourier transform of A (q): 

-in 



P(f) 



d 3 a 

e iqr A (a) — 



From the last three equations we can easily calculate P (r) (we use the nonrelativistic value W ~ 2m) : 

1 



1 d d 

P (f) = — (VlViPz) K*cr„^i) — — 

v A axi ox r , 



(16) 



(17) 



(18) 



1 



4r {{ax. a 2 ) - 3 (ai.n) (a 2 .n)) (1 + m A r) 



(ai.n) (a 2 .n) 

r 

In the last formulae we have omitted the spinors (fi t (p'^... but have labelled the cr's by the indices showing to which 
fermion they refer 

Following the discussion of the Introduction we shall consider only fermions which are heavy enough to bind. The 
radius of the bound states is of order a which must be small enough to ensure the validity of the inequality. 



rriAa <C 1 



So the potential ( |T8|) would be of the form: 



P(f) 



1 e" 



[c?i.<72 — 3 (ai.n) ((T 2 .n) 



(19) 



(20) 



To compare our nonrelativistic approach with relativistic equations we shall use the technique developed by 
Krolikowski in handling the spin dependence of the wave function B For the Schrodinger equation 



— V 2 tp + P (r)i> = Etp 
m 



(21) 



we decompose -0 into spin-zero and spin-one parts: 



where the projection operators are: 



1 



1 



P = -(l-ai.a 2 ), Pi = -(3 + ai.a 2 ) 



(22) 



(23) 
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The following set of identities can be used: 

Po (c?i - S 2 ) = {Si - S 2 ) Pi, P (Si + S 2 ) = (Si + S 2 ) P = 0, 



P (aii.a 2 n — 0i n cr 2 j) = iPo^int {?u - ail) 
Po (<Jn.a 2n + ai n a 2 i) — —26 in P , 

Po (fli — CT 2 i) ((Tin + <7 2 n) — 2iPo€i n i (oit — (J 2 l) 

Po (fii - a 2i ) (crin - a 2n ) = 4P 5 in 

Po (fii — &2i) (Pl n V2t + CU&2n) 

2 [S n e (an — <J 2l ) — 5 ni (<tu — &2i) — Su (a i n — <r 2n )] Pi 



Acting by Po on both sides of eq. (|2lj) one sees that the tensor forces vanish in the singlet state, so that <f> obevs the 
free Schrodinger equation. On the other hand by applying the operator (ai — S 2 ) Pi = Pq (Si — a 2 ) to eq. ( |2l| ) wc 
can derive the equation for the triplet \ function: 

-~V 2 X= 1 ^- T -^(x-3n(x.n)) = Ex (24) 

Following Q we now split x(r) into "electric", "longitudinal" and "magnetic" parts: 

X(r) =Xe (f) + XL (r) + Xm (r) , (25) 
where each obeys the following expansion in spherical harmonics: 

£(0 =^nT 3 - m (ni)xr(r), (26) 

jm 

d (Y jm j 



-Esc 



j i" 



*-«--£(«><£) MM™ 



In eqns. (Eq) we introduce differention in the transverse direction 



d d 

= (5 in - niiij) - — , (27) 

emu. ouj 

so that 

Ui^— = 0. (28) 

On the other hand the gradient and the Laplacian are expressed in terms of d/dr, d/dn± as follows: 

— --— +n— (29) 
dr r dnx dr ' 

, d 2 a 2 2 9 1 d 2 
V 2 = 



dr 2 dr 2 r dr r 2 dfi 2 ^ 
Note that d 2 /dn 2 L represents the angular part of the Laplacian. Note also that: 
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d 



d _ d 
ft = 2, 



n± I n±. 



dn ) ' dn\ 



Yjm = -J (.7 + 1) Y j7 



(30) 



The quantum number j corresponds to the total rather than the orbital angular momentum. Indeed let us define the 
total angular momentum J = L + S acting on the vector X of eq. (j2?3) as: 



{LXk)i = -ie ipq n p ^, (S X k) 



-itiktXl- 



(31) 



It is easy to see then that the operator L when acting on X e, XL and \m can be commuted through the vectors 
n, d/dn± and (n x d/dn±) and that the contribution of S in J = L + S cancels out so that J acts exactly in the 
same way as L on Yj m . For example, taking (\ e (r)) k = Xek (r), we get 



(JXek)i 



If we now repeat this operation: 



^ xl™ ( r ) l-itipqnp-J^- (n k Y jm ) - itiuniY ]r \ = 

■irn K ' Q ) 



(32) 



EX*™ W 



dY TI 



J 2 Xek = Ji ( J Xek), = J2 M L 2 Y Jm = J (j + 1) X ek, 



(33) 



which proves that j is actually the total momentum. The same proof can be given for Xl and xm ■ We are ready now 
to derive the equations for the radial wave functions. For that purpose we project eq. (24) on n and act on both sides 

Q 

of the equation by d/dn± and n x . Using the expansions (|25| ) and (|26| ) it is easy to obtain after some algebra 



the following equations: 

d 2 Xe 



2 d Xe J (j + 1) + 2 



dr 2 



r dr 



2xl fn 

Xe t - H FliXe + m£: Xe = 0; 



(34) 



dr 2 



dr 



-XL 



2i(i + i). 



c? 2 xm , 2dxAf i (i + 1) . 



dr 2 



dr 



-XM 



m 
2^1 



mxL 
2-Kv 2 A r z 



mExi 



Xm + mExM = 0. 



where we have omitted the indices jm of the functions xi"\ ... in these equations. 



We see that the functions Xe and xl are coupled in the first two equations ( |34| ) whereas the function xm is separated 
from XejXi- This is related to the fact that X M corresponds to the value of the orbital momentum £ = j while Xe 
and xl are the mixtures of the states x + with £ = j + 1 and x~ with £ = j — 1. These mixtures are [|J: 



2.7 + r 



XL = vW+1) 
Passing in eqns. ( |34| ) from Xe and xl to x^ we obtain: 

d 2 X " , 2d x - i(i-i) 



/2j + l X ' 
.7 + 1 _ , 

-x -+ 



2.7 + 1' 



2j + r 



dr 2 



dr 



mEx 



2irv 2 , L r 3 



3-1 . 
2.7 + 1' 



3y^7(7+l) , 
2 j + 1 ' 



0. 



d 2 X + { 2d X + {] + !)(] + 2) 
dr 2 r dr r 2 



■ mEx^ 



2irv 2 A r 3 



2j + r 



2j + l 



o. 



(35) 



(36) 
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We see from (^) that j (j — 1) in the first equation and (j + 1) (j + 2) in the second correspond to i (£ + 1) for 
I = j =F 1. Note that in each part of Eq. (|36| ) the potentials are attractive and singular as ~ 1/r 3 at r — > 0. Note also 
that xm in eqn ([54]) does not have this singularity. We show in the next Section that this singularity is not present at 
all in a fully-relativistic treatment of pseudoscalar exchange: in the Appendix we show that the apparent singularity 
arises only when a non-relativistic approximation is made. So Eq. ( |3^ ) is valid provided that r is not too small; more 
precisely that r > a, where the cutoff a is given by Eq. (|l^). 



III. THE RELATIVISTIC EQUATIONS 



To obtain a relativistic treatment of the problem and in particular, to understand the smoothing out of the singularity 
1 /r 3 we begin with the two-body Dirac equation (the Breit equation) 



£-7o 1} (l [1) V + mi 



(2) 

7o 



-i^p + m 2 ) - ^ (r) ^ = 



(37) 



describing a system of two spin-1/2 particles of masses mi and mi in the centre of mass frame, interacting with each 
other through a potential of the form 



Vint = V s (r) + V p (r) + V v (r) 

1) (2)' 

x 7o 



V s 



v„ 



(38) 



where S (r),P (r), and V (r) result from the exchange of scalar, pseudoscalar and vector particles. We use the Dirac- 
Pauli representation for the 7-matrices: 



7o=/3 



1 

-1 



7 



a 

-a 



75 = «7o7i7273 



(39) 



In this paper we shall be interested mainly in S and P exchange. Of course, a static potential V\ n t which depends only 
on the variable r = \f\ — f^l is not a relativistically invariant quantity. Furthermore the concept of a potential can 
only be an approximation at best to a proper field-theoretic description of a system. These are well-known difficulties 
of the generalisation of the Dirac equation involving a potential to of a two-body system which we do not wish to 
pursue here. We shall adopt the viewpoint that Eq. ( |37| ) provides a starting point for a calculation of the relativistic 
corrections to the non-relativistic two-fermion system in the same way that the treatment of the Coulomb potential 
in the Dirac equation leads to relativistic corrections in the spectrum of the hydrogen atom of higher order in a 
compared with the Schrodinger solution. An analysis of the constraints imposed by the relativistic invariance of a 
two-body fermionic system shows that provided the coordinates are chosen appropriately (which in practice means 
using the centre of mass frame for the equal mass case of mi = 1712), then a potential Vi n t (r) can be used in a 
relativistic two-body equation without violating relativistic invariance. Further details can be found in the papers 
by Mourad and Sazdjian || and Crater and Long ||, which also discuss the possible forms for the relativistic wave 
functions which result from this analysis. 



We now follow Tsibidis' 0] treatment of the Breit equation, which itself is based on the work of Krolikowski The 
spinor ip (r) is a sixteen component wave function which can be represented as a 4 x 4 matrix 



'0 '0 



ip- + ip — 



(40) 



where the indices +, — refer to the eigenvalues of 7q and 7Q 2 "*. Equation ( |37| ) is reduced to a set of 4 equations for 

"ip++, "0H Following the technique of ref. [3], which has been already used to analyse the non-relativistic case, we 

introduce the components 



(41) 



G 




where 4> and y correspond to the spin zero states while </> and \ correspond to S = 1 states (compare with eq. ( [22] ) . 
From Eqs. @), (|o|) and we now derive the following set of equations using the identities above for the projection 
operators Pq and Pi: 



\{E- 


s 


-p 


-V)<f>°- 


- ml + m2 (j) 


— ip.cj) 


= 


\{E- 


s 


+ p 


-V)4>- 


mi + m2 ^0 


= 




§(£ + 


s 


+ p 


-v)x°- 


m 1 -m 2 

2 A 


- W-X 


= 


+ 


s 


-p 


-v) x - 


mi—rri2 T ,0 
2 A 


= 






s 


-p 


-v)x- 


rai+ni2 -r?0 
2 A 


+ *PX° 


= 


\{E- 


s 


+ p 


-V)?- 


mi+ni! 
2 A 


+ ipx 


<f° = 




s 


+ p 


-V)$- 


mi-rn-2 ^0 


+ iptfr 


= 




s 


-p 


-V)$- 




+ ipx 


x° = o 



Each of the vector functions 4>, x, <j) , y°, entering this equations, can be split into "electric", "longitudinal" and 
"magnetic" parts according to Eqs. (|25|) and (|26| ) and expanded into spherical harmonics. We then obtain for each 
value of the total angular momentum 12 radial wave functions for the vector components : <f>> e (r) , 4> J L (r) , 4> J M (r) , x{ ( r ) 
etc, together with the four radial wave functions for the scalar components: OJ (r) , <j>> (r) , x°^ i r ) > X* i r )- We write 
down the sixteen equations for these functions according to the following classification. 

It is easy to see that (j),(jP , (/> e , 0° , (/>l , (/>° , xm , Xm have parity P = rj (— 1) J where r\ is the intrinsic parity (+1 for 
two fermions and —1 for fermion-antifcrmion system) and j is the total momentum (we omit index j at the radial 
components of the wave functions), the remaining 8 functions namely Xj X > Xe> Xe> XL, Xl> 4>Mi <^m' have P= —r\ {— 1) J 
and therefore the equations for the first 8 functions do not mix with the equations for the latter 8 components. The 
former case may be called a pseudoscalar meson trajectory (PMT) and its spectroscopic signature is 1 jj or 3 jj (the 
orbital momentum I = j) while the latter may be called a vector meson trajectory (VMT): it has the spectroscopic 
signature 3 (j — 1) ■ or 3 (j + 1) - g. A fermion-antifermion system which conserves charge conservation C will have 8 
independent components so the 16 spinor components will need to reduce to 8 dynamical equations and 8 constraint 
equations, or 4 dynamical equations and 4 constraint equations for PMT and the same again for VMT. We demonstrate 
this below. 

Thus we obtain after some algebra: 
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(i) pmt 3 j j ,p= v (-iy 



1 , mi + mo n 

-(E-S + P-V)^.- ' 2 V = 

I( E - S -P-^ )x „-!=i±^ x O M = 

i (£ - s + p - v) A - ^l XM + iM±n<i + (A + i) « _ o 

i (£ + 5 + P - V) L - ^^^< </ - ' < = n 

i(i, + 5-P-V)A-^,-(l + iU = 





2 


1 


mi 




m 2 




2 




mi 




m 2 




2 


< 


mi 




m 2 



(43) 



(ii) VMT 3 (i±l) j ,P=-7 ? (-l) 



(44) 



1 /t-, ^ n Trt n mi — m 2 ( d 2\ 1 
-(E + S + P-V) X ° -_ x -^_ + -j Xe -_ Xi = 

^ + S-P-^-^pV = o 

i(,-S-P-^e-^X° + f =0 

i( S -5-p-y) Xi -^^xi-^^x = o 

i(^, + P-^) x o_- + - Xi _^ + i)^ = 
i(^ + 5 -P-F)^_!^ 0M + i^ x o_^ + i )x o =o 

i( J E + 5 + P-^)^-^^^ M =0 

The non-relativistic analysis of the previous section for pseudoscalar exchange refers to the spin S = 1 states 

connected to the "large" components of the relativistic wave function -0++- This implies that \e, XL of Eq. (|||) are 
the non-relativistic analogues of the same functions entering Eq. (|44|) . Remember also that the 1 /r 3 singularity of 
Eq. ( |34| ) only involves the functions Xe and \l, so it is sufficient for our purposes to concentrate on the following set 
of equations for the case mi = m 2 = m: 
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(E + S + P - V) X " - 



d 2 
dr r 



X e - -XL = 



-(E-S-P-V) X e-mx° e + -£ r =0 
1 



-{E-S-P-V) X L-mxl- 



(45) 



{E - S + P - V) xl ~ mXe + 



(E-S + P-V)x° L -mxL- 



d 1 

dr r 



(E + S-P-V) cjP M 



i(i + i) v o 

Xe 



d l 

dr r 



X° L =0 



First of all we should try to show that the non-relativistic limit of ( f45|) gives Eq. (|3J) for x e and xl- This is not 
a trivial exercise but we can proceed as follows. Add and subtract the second and the fourth equation of (|44|) and 
introduce xt — \ (Xe ± xT)- Then we can get the equation (keeping only P ^ 0): 



(£7-2m)x+ + ^(l 



P 



E + 2m 



P 



E + 2m 



P 2 



E + 2m 



X t 







From ( [45| ) and (46) we can obtain for xt 



[XL+Xl) 



(E - 2m 



P 2 



E + 2m 



xt = o 



(46) 



(47) 



In the nonrelativistic approximation x°and <p° M are readily expressed through x+and xt, (this is correct only for the 
non-relativistic limit). We get: 



t>°M 



77! 
1 

77! 



(dxi 



\ dr r 
3 (J + 1) 



-Xl 

d 1 

dr r 



1 - 



P \ P' 



Xl 



4m J 4m 
P 
4m 



X,' 



(48) 



1 



EL 

4m 



xl 



where P' = In these equations we neglect all non-relativistic corrections and keep only linear terms in P. Then 



substituting ( ^8|) into (|i^) and ( f47|) we see that only derivatives of P remain in the final equations for x+and Xl- 
Eventually we obtain the non-relativistic limit: 



m{E — 2m) xt + 



. , d 2 xt , 2rf Xe + ]{] + !) + 2 



dr 2 



dr 



>X L 



P" 



P' 



4m 8mr 

, n „ s , rf 2 xl 2rfxl + 4. 2? (i + i) p" , 

m(P-2m)xU^ + l4 L - lil ^^-^^X+- i ^xJ = 



0, 



(49) 



dr 2 r dr 
For the exchange of a pseudoscalar particle: 



771 1 

P=—-e~ mAr 
ttv < , r 



(50) 



So we see that the terms proportional to P'/r and P"in Eq. ( [49| ) give the 1/r 3 singularity of eq. (f34|). We show in 
the Appendix that the relativistic equations of Eq. (Q) do not give a 1/r 3 singularity. Thus this singularity is just 
an artefact of the non-relativistic approximation. 
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IV. CONCLUSIONS 



We conclude that heavy fermions / of mass larger than 800 GeV are required if they are to bind with / by Higgs or 
pseudoscalar exchange. If fermions of this mass exist, a whole new region of physics will be associated with this high 
mass scale: the fermions themselves and their decay modes, and the / / bound states, their decays and systematics. 
The details of this new physics, which involves binding by an unfamiliar mechanism, cannot be predicted in advance. 

We have also demonstrated that the 1/r 3 singularity seemingly present in pseudoscalar exchange will vanish if a 
relativistic framework is adopted for the calculation. Pseudoscalar exchange in relativistic quantum field theory arises 
from renormalisable theories and hence cannot lead to collapse. 

We would like to thank G. Tsibidis for his help with the Breit-Krolikowski equations. This work has been supported 
in part by INTAS grant 93-283 and INTAS-REBR 943986. One of us (AA) thanks PPARC and the University of 
Sussex for support during his stay at Sussex in 1997. 



V. APPENDIX 



We now look at the relativistic case to show that there is no 1/r 3 singularity at the origin arising from pseudoscalar 
exchange. It is not easy to derive relativistic equations for the two functions % e and xl of the non-relativistic ap- 
proximation (E^). It is easier to exclude x&i Xe, XL, Xl from Eq. (fi"5|) and then obtain a pair of equations for x°i 



To this end we have from the second and fourth equations of (Eq) (keeping P ^ and S ^ 0): 



Xe 



Am 



Am 2 -(E- S) z + P 2 
4m 



4m 2 - (E - S) + P 2 
and from the third and the fifth equations: 



1 

2m 

1 
2m 



{E-S + P) 







dr 


r 






r 


dr 



(51) 



-4m 



XL 



Xl 



4m 2 - (E - Sy + P 2 

—4m 

4m 2 - (E - S) 2 + P 2 



2m v 



2m 



(E-S-P) 



d 
dr 



d_ 1 

dr r 



(52) 



Substituting these expressions into the first and the last of equations ( [45]) we get two second order equations for x°and 



d 2 X n , 2d X ° J(j + 1) o 
dr 2 r dr r 2 



d 
dr 



In 



4m 2 -(E- Sy + P 2 



4m 



E-S + P 
d 



Am 2 

dx^_ 
dr 



(e - sy + p 2 



S + P 



S + P 



X" 



E-S + Pdr 



d 2 Ki 
dr 2 
_d_ 
dr 



2d(jP 



M 



In (4m 2 -(E-Sy + P 2 
j (3 + 1 



t>°M 



In- 



r dr 
4m 2 — 



(e - sy + p' 2 



In [Am 2 -{E-Sf + P 



E-S-P 

2 . ,v2 



Am 2 

— + - 

dr r 

1 



(E - Sf + P 2 
Am 



E + S-P 



E-S 
d 



E - S-Pdr 



x° = o 



(53) 



It is now clear what a delicate problem it is to go to the nonrelativistic limit for the pseudoscalar case. For the scalar 



interaction, S should be retained only in 4m 2 — (E — Sy 



-Am (e — S) where e = E — 2m in the first line of Eq. 
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(|53|). This results in the usual Schrodinger equation with potential energy S. For S = 0, P 7^ we should keep the 
apparently small term P/2m even in the limit m — > 00 term since the coupling constant in P is proportional to m 2 . 
We do not however keep P 2 /4m 2 terms in order to get the nonrelativistic equations (36). The most important thing 
which we learn from ( p3|) is that when r — ► there are only singularities ~ 1/f 2 , and therefore there is no collapse. 
A more detailed analysis of the relativistic equation is outside the scope of this paper: Further details can be found 
in the papers by Mourad and Sazdjian Q and Crater and Long ||. 
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